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Abstract. Wc address the problem of obtaining the exact reduced dynamics of the 
spin-half (qubit) immersed within the bosonic bath (enviroment). An exact solution of 
the Schrodinger equation with the paradigmatic spin-boson Hamiltonian is obtained. 
We believe that this result is a major step ahead and may ultimately contribute to 
the complete resolution of the problem in question. We also construct the constant 
of motion for the spin-boson system. In contrast to the standard techniques available 
within the framework of the open quantum systems theory, our analysis is based on 
the theory of block operator matrices. 
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1. Introduction 

The Hamiltonian of the paradigmatic spin-boson (SB) model is specified as [1-4] 

H SB = H s ® I B + I s ® H B + H int , (1) 

where 

POO 

Hs = (/3(T Z + oca x ) and Hb = / dujh(uo) a)(uo)a(uj), (2) 

Jo 

are the Hamiltonian of the spin-half (qubit) and the bosonic field (environment), 
respectively. The interaction between the systems has the following form 

poo 

Hi llt = (? z ® / dtu (g(u)*a(u) + g{u)a\uj)) = a z ®Y. (3) 
Jo 

Is and Ib are the identity operators in corresponding Hilbert spaces of the qubit and 
the environment, respectively. 

In the above description, a z and a x are the standard Pauli matrices. The bosonic 
creation a'(cj) and annihilation a(cj) operators obey the canonical commutation relation: 
[a(u) , at (r])] = S(cu — t])Ib, for w,?| e [0, oo). The functions h, g G L 2 [0, oo] model the 
energy of the free bosons and the coupling the bosons with the qubit, respectively. The 
constants a and /3 are assumed to be real and non negative numbers. Furthermore, (3 
represents the energy gap between the eigenstates |0) and |1) of a z , while a is responsible 
for the tunneling between these states. The Hamiltonian (1) acts on the total Hilbert 
space Titot = C 2 ® Tb, where J-b := ^ r (-^ 2 [0, oo]) is the bosonic Fock space [5]. 

It is worth mentioning that more often we encounter situations in which there is 
a countable number (finite, in particular) of bosons (see e.g. [6-9]). In such cases we 
define the SB model via the following Hamiltonian 

H SB = (/3cr z + aa x ) <g> I B + I s ® h k a[a k + a s ®y] (g* k a k + g k a ] ) , (4) 

k k 

where the creation and annihilation operators a k , a k satisfy [a k , a\\ = 5 k [. Formally, it 
is possible to obtain (4) from (1) by setting 

x(u) = x k 5(cu — cu k ), where x = h,g. (5) 

k 

Therefore, we can treat both cases simultaneously. Although generalizations of the SB 
model (e.g. asymmetric coupling [10]) are also under intensive investigation, we will not 
focus on them in this paper. 

The problem of a small quantum system coupled to the external degrees of 
freedom plays an important role in various fields of modern quantum physics. The 
SB model provides a simple mathematical description of such coupling in the case of 
two-level quantum systems. For instance, an interaction between two-level atoms and 
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the electromagnetic radiation can be modeled via the SB Hamiltonian [11]. For this 
reason the SB model is of great importance to the modern quantum optics. There are 
various physical problems (e.g., decoherence [12-15], geometric phase [16]) related to the 
properties of the model in question, which has already been addressed and intensively 
discussed. Nonetheless, an exact solution of the Schrodinger equation: 

id t \%) = H SB \%) with |* ) = |*}, (6) 

is still missing for both a ^ and f3 ^ 0. Several approximation methods [17] have 
been developed in past fifty years to manage this problem. Models obtained from the 
SB Hamiltonian under mentioned approximations are well-established and in most cases 
they are exactly solvable. The famous Jaynes-Cummings model [18] can serve as an 
example. Formally, one can always express the solution of (6) as |\I/ t ) = Ut|^), where 
Ut := exp(— iHsB^) is the time evolution operator (Stone theorem [19]). Needless to 
say, such a form of the solution is useless for practical purposes. 

There is at least one important reason for which a manageable form of the time 
evolution operator XJ t is worth seeking for. Namely, it allows to construct the exact 
reduced time evolution of the spin immersed within the bosonic bath, the so-called 
reduced dynamics [20]: 

p t = Tr B (U t p ®w B Ui). (7) 

Above, the state cj b is an initial state of the bosonic bath. Tr B denotes the partial trace, 
i.e., Tr B (M ® X) = MTrX, where Tr refers to the usual trace on Tb- For the sake of 
simplicity we have assumed that the initial state of the composite system p int is the 
tensor product of the states po and w B . In other words, no initial correlations between 
the systems are present [21-24]. 

In general, the formula (7) is far less useful, than its theoretical simplicity might 
indicate. Indeed, to trace out the state U t pi n tUj over the bosonic degrees of freedom, 
one needs to i) calculate U t and ii) apply the result to the initial state p int . Herein we 
will cover the first step and we will investigate the ability to accomplish the second one. 

In order to write the time evolution operator U t in a computationally accessible 
form, the diagonalization of its generator Hs B or an appropriate factorization [25] is 
required. It can be found (see e.g., [26-28]) that the problem of diagonalization on 
the Hilbert space C 2 ® can be mapped to the problem of resolving the Riccati 
equation [29]. This new approach was recently successfully applied to the time- 
dependent spin-spins model [30]. As a result, the exact reduced dynamics of the qubit 
in contact with a spin environment and in the presence of a precessing magnetic field 
has been obtained. It is interesting, therefore, to apply this approach to the SB model 
as well. This paper has been devoted to accomplish this purpose. Although, an explicit 
form of the Riccati equation has already been derived [31], the solution has not been 
provided yet. In this paper we derive an exact solution of this equation assuming (3 = 0. 
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2. The block operator matrix representation and the Riccati equation 



We begin by reviewing some basic facts concerning a connection between the theory of 
block operator matrices [32] and the SB model. First, the Hamiltonian (1) admits the 
block operator matrix representation [31,33]: 



H 



SB 



H B + V + (3 
a 



Hi 



Oi 

- V 





H + 


a 




a 


H 



(8) 



with respect to the direct sum decomposition H tot = J-"b© J-b of Utot- The entries a and 
(3 of the operator matrix (8) are understood as ale and /31b, respectively. Henceforward, 
we use the same abbreviation for any complex number. 

The Riccati operator equation associated with the matrix (8) reads [31] 



aX 2 + XH + - E_X 



a 



0. 



(9) 



where X is an unknown operator, acting on J 7 ^, which needs to be determined. The 
solution of this equation, if it exists, can be used to diagonalize the Hamiltonian (8). 
To be more specific, if X solves (9) the following equality holds true 



S HsbS 



E + + aX 
H 







a 



where S 



1 

X 



1 



(10) 



By means of this decomposition we can write \J t in an explicit matrix form: 

U, = Sdiag[e- i ( H ++^,e- i ( H -- frft >]S- 1 . 



Note, the last formula reduces the problem of finding the solution of the Schrodinger 
equation (6) to the problem of resolving the Riccati equation (9). It is well-established 
that the reduced dynamics (7) can easily be obtained when a = [5]. In this case 
no additional assumptions on [3 are needed, which should not be surprised since the 
matrix (8) is already in a diagonal form (X = 0). Moreover, if a = the qubit does 
not exchange the energy with the bosonic field because [Hs <8> Ib, Hsb] = 0. Therefore, 
the only exactly solvable case, which is known at the present time, represents rather 
extreme physical situation. 

In the next section we will derive an exact solution of the RE (9) assuming (3 = 0; 
nevertheless, we not impose any restrictions on a. This is exactly the opposite situation 
to the one we have discussed above. At this point, the natural question can be addressed: 
what about the case, when both a and /3 are not equal to zero? Unfortunately, the 
answer is still to be found. In fact, usually the SB model is defined only for /3 = 0. At 
first, it might seem that the complexity of the problem is the same both for (3 = and 
(3 ^ 0. Although, this is indeed true when a = 0, no argument proving this conjecture 
for a 7^ has been given so far. We will return to this matter at the end of the next 
section. 
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3. Solution of the Riccati equation 

3.1. Single boson case 

To understand the idea of our approach better let us first consider the case where there 
is only one boson in the bath [34,35]. Then, the Hamiltonian of the SB model can be 
written by using the block operator matrix nomenclature as (/3 = 0) 



H 



SL! 



H a 

a H + 



with H ± = <jja}a ± (g*a + gcv) 



The operators H ± can be expressed in a more compact form, that is 



(12) 



H = uj~D fa) aD_f — E and H + = tuD^fa^oD y — E, 

where / = g/u and E = \g\ 2 /oj. The displacement operator D/ := exp(/*a 
the following, easy to prove, properties 



(13) 
fa)) has 



i) D_ / = D}, ii) D/D.^Ib and hi) D f D g = e & ^D f+g . (14) 

Q stands for the imaginary part of the complex number fg*. The relations (13) can 
be proven by using equality D^aD_j = a — f, which follows from the Baker- Campbell- 
Hausdorff formula [36, 37]. For the sake of simplicity and without essential loss of 
generality we rescale the Hamiltonian (12) so that Hsb - > Hsb + E. This is nothing 
but a rescaling of the reference point of the total. 

After this procedure the Hamiltonian (12) takes the form 



H 



SB 



a LoD_fa>aDf 
while the corresponding Riccati equation reads 

aX 2 + X (uD f a^aD_ f ) - (cuD_ f a^aD f ) X - a = 0. 
To solve this equation, let us first define an operator in a way that 



(15) 



(16) 



P v := exp(i^a^a), tp £ [0,2n). 



It is not difficult to see that 



;i7) 



pt 



p p 



Ib and iii) P^P 



Moreover, from the Baker-Campbell-Hausdorff formula we also have P^aP. 
which ultimately leads to 



(18) 



P ¥ ,D / P_ V 



(19) 
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In what follows, we will prove that P„- solves the Riccati equation (16). First, let 
us note that P„. is a function of the number operator a^a, thus [P„., a^a] = 0. In view 
of (19) we obtain P 7r DjP_ 7r = D_j, hence 

P n (D f JaD_ f ) = (D_f<JdD f ) P ff . (20) 
By writing P n in terms of the eigenstates \n) of a^a we obtain 

P w = £e*"»<n| = ^(-1)»R (21) 

ngN ngN 

where we used the well-known mathematical fact that a)a\n) = n\n), for n G N. Finally, 
from (21) we conclude that P„. is an involution, i.e., P^ = Ib, which together with (20) 
leads to 

aP 2 n + P n (uD f a j aD_ f ) - (wD./a+oD/) P n - a = 0. (22) 

Note, P n transforms the creation and annihilation a operators into —a} and —a, 
respectively. In other words, Pjr can be interpreted as the bosonic parity operator [38]. 
Moreover, P„. does not depend on the parameter a; in particular, P n remains a nontrivial 
(X 7^ 0) solution of the Riccati equation (16) even when a = (Sylvester equation). 

Now, by means of the parity operator P„., we can derive an accessible form of the 
time evolution operator U*. According to (10) and (11) we have 

u, = i 

2 

where the quantities U±(i) and V±(£) read as follows 



U + (t) V + (t)P w 
V_(t)P^ U_(t) 



(23) 



XJJt) = e ~ i ( U ±+ aP ^) t + e -i(H±-aP^)^ YJt) = e -i(H±+aP,r)t _ e -i(H±-aP^)t_ ^4) 

For a = the formula (23) simplifies to the well-known result [5], which can be obtained 
independently, without solving the Riccati equation. 

It is instructive to see how the bosonic parity operator P^ can also be used to 
construct the constant of motion for the SB model. For this purpose let us take 
J-tt : = &x ® Pn', then [J n , H SB ] = 0, thus from the Heisenberg equations of motion 
follows J n = 0, which means that J n does not vary with time. Since P w is an involution, 
i.e., P^ = Ib thus 3 n is an involution as well. Therefore, J n can be seen as the parity 
operator of the total system. In conclusion, the total parity is conserved when (3 = 0. 

For (3 7^ the parity symmetry of the total system is broken and the Riccati 
equation (16) cannot be solved by applying a similar method to the one we have 
used above in the case of /3 = 0. From mathematical point of view, the problem 
arises because the diagonal entries He ± V ± (3 are no longer related by an unitary 
transformation. Indeed, if the converse would be true, there would then exist an 
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unitary operator W such that (H B + V + j3) W = H B — V — (3. Thereby the spectra 
a (H B ± V ± f3) = cr(H B ± V) U {±/3} would be the same, which clearly is impossible 
unless (3 = 0. As a result, for q ^ one can expect that the mathematical complexity 
of the SB model is different within the regimes /3 = and (3^0. 

3.2. Generalization 

The results of the preceding subsection can be generalized to the case where there is 
more that one boson in the bath. In order to achieve this objective one needs to redefine 
the displacement operator D/ in the following way 

D f ->■ exp (A - A ] ) , where A = V ^a k . (25) 
Then, the solution of the Riccati equation reads 

P = exp(i7T ^ 4a fc ) = (%) P^fc, where P^ = exp^a^a*). (26) 
k k 

4. Remarks and Summary 

In this article, we have solved the Riccati operator equation associated with the 
Hamiltonian of the paradigmatic spin-boson model. Next, in terms of the solution 
we have derived an explicit matrix form of the time evolution operator of the total 
system. This, in particular, allows us to solve the Schrodinger equation (6). We wish to 
emphasize that in order to obtain the reduced dynamics (7) one more step is required. 
Namely, the terms like 

Tr(e- i[H±±aP - ] W i[H±±aIV1 ') (27) 

need to be determined. Of course, one can always evaluate the quantities given above 
by using e.g., perturbation theory. However, the true challenge is to establish this 
goal without approximations. It seems that the simplest way to do so is to solve the 
eigenvalue problem (H ± ± aP n )\ip) = The ability to solve this eigenproblem 

separates us from deriving the exact reduced dynamics of the qubit immersed within 
the bosonic bath. We stress that for a ^ the problem is nontrivial since the qubit 
exchange the energy with its environment. Moreover, an impact on the mathematical 
complexity of the model, has not only a transfer of the energy between the systems but 
also the energy split between the states |0), |1). 

Interestingly, the Riccati equation is a second order operator equation, thus one 
can expect that its solution involves a square root. In particular, nothing indicates 
that the solution should be linear as it is in our case. Therefore, we not only solved 
the Riccati equation (9) but also linearized the solution. At this point, a worthwhile 
question can be posed: is it a coincidence that the linear operator happens to solve 
a nonlinear equation? Perhaps, it is a manifestation of some additional structure in 
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the model. Historically, a similar situation took place when Dirac solved the problem 
with a negative probability by introducing His famous equation [39]. By linearizing the 
Hamiltonian of the relativistic electron, Dirac not only predicted an existence of the 
antiparticles but also explained the origin of the additional degree of freedom of the 
electron. 
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